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Abstract. To investigate the distribution of dust in Sb and 
Sc galaxies we have analyzed near-infrared and optical surface 
photometry for an unbiased sample of 37 galaxies. Since light 
7— H in the A'-band is very little affected by extinction by dust, the 
B — K colour is a good indicator of the amount of extinction, 
and using the colour-inclination relation we can statistically 
determine the extinction for an average Sb/Sc galaxy. 

We find in general a considerable amount of extinction in 
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spiral galaxies in the central regions, all the way out to their 



effective radii. In the outer parts, at Djf ,21 , or at 3 times the 
typical exponential scale lengths of the stellar distribution , we 
OO find a maximum optical depth of 0.5 in B for a face-on galaxy. 
^— 5 If we impose the condition that the dust is distributed in the 
same way as the stars, this upper limit would go down to 0.1. 

O 
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1. Introduction 

C^The analysis of surface photometry in B of large samples has 
^ shown that Sb and Sc galaxies suffer from much more extinc- 
tion than was previously thought (Valentijn 1990). In the cen- 
"jrj^ tres there is so much extinction that generally the dust optical 
depth in the visible is larger than unity. There are claims that 
these spiral galaxies are optically thick even in the outer regions 
(Gonzalez-Serrano & Valentijn 1991; Burstein et al. 1991; Bot- 
tinelli et al. 1995; Valentijn 1994, hereinafter V94) although 
these are disputed by several authors (White & Keel 1992; 
Huizinga & van Albada 1992). 

In this letter we present the results of a direct way to mea- 
sure the amount of extinction using radial surface brightness 
profiles in B and K for a sample of galaxies. Since the ex- 
tinction in K is generally more than an order of magnitude 
lower than in B, a change in B — K , in the absence of a stellar 
population gradient, indicates only a change in the amount of 
extinction. Here, we separate the effects of extinction from stel- 
lar populations by using the fact that extinction is inclination- 
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dependent, while radial stellar population gradients are not. 
Thus, by taking a sample of galaxies uniformly distributed in 
inclination, we can determine statistically the amount of ex- 
tinction at various places in the galaxy. This cannot be done 
for optical-optical colours because they suffer from the fact 
that not only the colour of the underlying stellar population 
is unknown, but also the extinction law, since dust and stars 
are mixed, so that the effective extinction law depends on the 
quantity and the distribution of the dust (see e.g. Evans 1994; 
Jansen et al. 1994). We present here a study of 37 galaxies, for 
which the data were presented in Peletier et al. (1994, Paper I). 
The data allow analysis out to the 21 st A'-band isophote, which 
is close to the 25 th 5-band isophote. We analyse in several dif- 
ferent ways the results in the centre, at an effective radius (the 
radius inside which half the light is contained), and at Djf,2i, 
and discuss the consequences for the optical depth of galaxies 
at these locations. 

2. Optical-infrared colour distributions 

The sample analyzed here consists of galaxies of morphologi- 
cal type 3-6.4, corresponding to Sb and Sc, with relatively high 
central Blue surface brightness (20.5-21.5), selected on the ba- 
sis of the ESO-LV catalogue (Lauberts & Valentijn 1989), i.e. 
normal late-type spiral galaxies of the types that contain large 
amounts of extinction (Valentijn 1990). The sample is diame- 
ter limited and uniformly distributed in orientation on the sky. 
A detailed description is given in Paper I. 

The advantage of having extinction-free K-band photometry 
is that it allows us to define a number of radii in each galaxy 
that depend only on the true stellar distribution, unaffected 
by extinction. For example, when analyzing the colours at an 
effective radius in K , we do not have to worry how much the 
measurement of the effective radius itself at various radii is 
affected by the extinction. We have defined 3 radii at which 
we study the dependence of B — K as a function of axis ratio. 
The colours are given in Fig. la. On the left we plot the ex- 
trapolated disk colour (io,B - t^o,K, in the middle the colour as 
measured at the effective radius in K and, in the right panel, 
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21. Fig. 1 shows the following effects: 
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Fig. 1. Azimuthally averaged B — K colour at various positions 
in the galaxy. On the left is shown the extrapolated central disk 
colour, in the middle the colour at Y e j / ( A', and on the right at K=21 
mag arcsec - 2 . Galaxies for which no redshift is available, and for 
which as a result Mj( could not be determined, are indicated with 
asterisks. Filled symbols are given to galaxies brighter than Mj^ = 
—22, while fainter galaxies are indicated with open symbols. The lines 
show the colour-dependence as a function of inclination for simple 
models with central face-on optical depths r^ ac€ of 0, 0.5, 1, and 5. 



i) We find that these galaxies are much redder in the centre 
than in the outer parts. A typical difference in B — K between 
the centre and Djf,2i is one magnitude, in some cases more. 
This is much more than seen in, e.g., elliptical galaxies. 

ii) In the centre, the spread in B — K for face-on galaxies 
is much smaller than for edge-on galaxies. Some edge-ons are 
much redder than all face-ons, while others are bluer. Fig. 1 
shows that the edge-on galaxies that are bluer than all face-on 
systems are either fainter than Mjf = -22 or have unknown 
redshifts (and so could also be fainter than Mjf = -22). Not 
all faint systems are edge-on, but the bluest are. We think this 
is because these galaxies have been misclassified. Edge-on sys- 
tems are much harder to classify usually than face-ons: small 
edge-on Sb's on a plate do not look very different from edge- 
on Sd's. We argue that most of the blue edge-on galaxies here 
have a type later than 6.5 and have been misclassified in the 
ESO-LV catalogue (Lauberts & Valentijn 1989). This justifies 
performing a second selection on the sample, using the ii'-band 
luminosity Mjf in order to obtain a more homogeneous sample 
for further inclination test studies. (As a matter of fact, one 
of our conclusions from the K-band photometry is that any 
sample defined on the basis of its K-band properties is much 
cleaner than a sample selected on the basis of visual passband 
photometry. [Paper I].) After having selected \Mk\ > 22 the 
situation is much clearer. We find that in the centre (filled dots) 



and at r e ff edge-on galaxies generally get redder strongly with 
inclination whereas at Djf,2i no trend of B — K as a function of 
inclination can be determined. In Fig. 1 we have also displayed 
the behaviour of a number of models. These are simple models 
(see Section 3) with a range in face-on extinction. In the central 
regions, the majority of the galaxies follow the models with a 
central face-on optical depth of 2 or more. In the outer parts, 
the transparent models seem to be preferred. One should how- 
ever note that even after applying the luminosity criterion the 
range in colour for edge-on galaxies is still much larger than 
for face-ons. This clearly shows that at a given inclination the 
amount of extinction varies from galaxy to galaxy. The limit 
at Mjf = -22 is not critical. The relations do not change very 
much after applying a different Mjf selection. 

We conclude from Fig. 1 that galaxies seem to behave as 
expected if they are opaque in the centre and at r e ff but are 
transparent at Djf ,21 • 



3. Analysis of the scalelength ratios 

In this section we use the slopes of the radial surface brightness 
profiles, or the scale lengths, to determine the amount and the 
distribution of the dust. We consider a galaxy that is represen- 
tative of the sample as a whole: the ratio of exponential scale 
lengths in B and K (Xb/<xk equals 1.3, when face-on, and 1.7 
when seen edge-on (Fig. 4 of Paper I). 

The fits to obtain these scale lengths have been made on the 
radially averaged surface brightness profiles (Paper I), between 
1 and 3 ii'-band scale lengths. Here, exactly the same radial 
range has been used for the fits in B and K . First we have to 
correct for the effect of radial abundance gradients. Two inde- 
pendent lines of reasoning are used here (see Paper I). First, for 
galaxies similar to the ones discussed here, but without much 
dust, i.e. SO or Sa galaxies, Balcells & Peletier (1994) find that 
the ratio of the scale lengths in B and I equals 1.04 ± 0.05. 
Using almost any stellar population model this gives a scale 
length ratio between B and K of 1.08 ± 0.10. The second ar- 
gument comes from abundance gradients measured from Hll re- 
gions. Literature data (Vila-Costas & Edmunds 1992; Zaritsky 
et al. 1994) indicate that the average gradient in normal nearby 
galaxies equals A(0/H) = —0.25 dex/(5 scale length), or - 
0.17 dex/(i<f scale length). This corresponds to a scale length 
ratio as/ajf = 1.17, using a simple single-age stellar popula- 
tion model. This corresponds more or less to the lower limit 
of our observations. We assume in the rest of this letter that 
stellar population gradients are responsible for a scale length 
ratio as/etK = 1.1, independent of axis ratio. So the dust will 
only have to explain as/etK = 1.2 for a face-on galaxy, and 1.6 
for an edge-on. De Jong (1995) claims that the stellar popu- 
lation gradients in these spirals are much larger than assumed 
here. His reasoning is that scattering by dust is so important 
that the amount of extinction needed to explain the gradient 
in a blue colour like B — V would imply a gradient in B — K 
which is much larger than observed. Because of the discrep- 
ancy with metallicity measurements from Hll regions he also 
needs considerable age gradients. His explanation may be valid 
for his data, consisting only of face-on galaxies, but with his 
inferred stellar population gradients one would not expect the 
good correlation with inclination that we see here. Given also 
the fact that there would be a discrepancy with SO's and Sa's 
(see above) we suggest that the explanation given here is more 
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realistic. 

To model the observed scale length ratios we first tried to 
use the relation of Valentijn (1991) between a b s , ce* and ad, 
the observed scale length, the stellar scale length, and the dust 
scale length. The model gives an analytic solution to a radial 
version of the equation of radiative transfer. It includes uniform 
distributions of dust and stars in separate layers. We found that 
the results depend quite strongly on the vertical distribution of 
the dust and stars, and since there is a large difference between 
uniform and exponential vertical distributions, we resorted to 
numerical modeling. Our exponential dusty galaxy model con- 
tains a uniform distribution of stars, exponential in the radial 
and vertical direction, and a uniform distribution of dust, also 
exponential in both directions, but with different scale lengths 
and heights. The ratio of the scale heights of dust and stars we 
call £. We project the models and fit exponentials to the major 
axis profiles in the same way as is done for the observations. 
One set of models has been made for face-on galaxies and one 
for galaxies at 80 degrees. 

The scale length ratios are plotted in Fig. 2 for various values 
of £ and ctd/ctK- We impose (driven by the observations) the 
condition that as/ajf in the face-on case lies between 1.1 and 
1.3, and in the edge-on case between 1.5 and 1.7. Furthermore 
T edge should be 5 times as large as r^ ace , reflecting the differ- 
ence in path length through the dust. Table 1 indicates which 
model galaxies fulfill the first criterion (column 3), which the 
second (column 4), and which ones fulfill the latter as well, i.e. 
are acceptable solutions (column 5). 
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We find solutions for any thickness, if ad/ a* > 2. The central 
face-on optical depth is around 1 or 1.5. For a model with 
«d/«A'=2, we would expect tb=0.2-0.3 at 3 scale lengths for 
a face-on galaxy, or tb = 1 if the galaxy were edge-on. It is 
interesting to see that no solutions are found with ced/ce* = 1. 

The current model is too simplistic, especially because it 
ignores scattering, and because of the fact that dust is generally 
distributed non-uniformly across the galaxy. 

The effects of scattering have been modeled by e.g. Witt 
et al. (1992) and Byun et al. (1994). The paper by Witt et 
al. shows that although it is difficult to calculate the effects 
of scattering one can describe it to a good approximation by 
modifying the extinction law, i.e. by decreasing the absorption 
in the optical relative to the IR (see also Rix & Rieke 1993). 
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Fig. 2. Scale length ratios between B and K as a function of central 
optical depth for disk models with various ratios of scale lengths and 
heights of dust and stars. 



Scattering is more efficient for face-on galaxies than for edge- 
on's since, in the first case, there is a larger probability that 
photons are scattered out of the plane towards the observer. 
Changing the extinction law by means of including the effects 
of scattering implies that more dust is needed to explain the 
behaviour in B (a factor 2 in optical depth at most). The fact 
that scattering is more efficient for face-ons also means that a 
model including scattering requires less extinction for edge-on 
galaxies relative to face-ons than a model without scattering. 

Let us now assume that the dust is distributed in clumps, as 
opposed to being uniform across the disk. We are interested in 
the case in which a substantial fraction of the disk is covered 
with dust clouds. If the fraction is much less, scale lengths are 
not affected much, contradictory to the current observations. A 
galaxy with a face-on cloud coverage of ~ 50% has an edge-on 
coverage of close to 100%. This means that the extinction is 
much more efficient for the edge-on than for the face-on since, 
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in the latter case, never more than 50% of the light is being 
extincted. This means that the edge-on optical depth needed 
to produce the observed scale lengths is less than 5 times the 
face-on optical depth. How much less depends on the geometry 
of the dust. 

When both effects are included, the models will need more 
face-on extinction, and a smaller ratio for face-on to edge-on. 
Taking as an example the solution for which T ed 3 e j T f ace = 
2.5, and multiplying the face-on optical depths of Table 1 by a 
factor 2, the model is in agreement with the observations for a 
central face-on optical depth in B between 2 and 3. 

4. Discussion 

In the previous section we have found that the only solutions 
that agree with the observed scale length ratios in B and K 
are solutions with large ratios between the scale lengths of dust 
and stars. Including absorption alone we do not manage to find 
a solution with equal scale length of dust and stars. 

However, including scattering and dumpiness in the analy- 
sis, it is possible that one will find acceptable solutions with 
ad/a* = l. These will have a central optical depth in B of about 
2. In the outer parts, at 3 stellar scale lengths, the face-on op- 
tical depth would be reduced to 0.1. 

What is the maximum amount of extinction at Djf,2i? Here 
the scale length model of section 3 is less useful than the di- 
rect observations of section 2. Fig. 1 shows that models with 
T l\ Ce ^ 0.5 have big problems fitting the colours, especially 
when turned edge-on. Since the data here do not exclude rl^ ce 
between 0.1 and 0.5 it might be attractive to make a two- 
component model of the dust (see Valentijn 1995). A one- 
component model with a large ratio for ced/ce* (e.g. 3) would 
either have too little extinction in the centre or too much in 
the outer parts. Parameters of an acceptable two-component 
model are given in Table 2. 
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An advantage of having these two-component models is that 
diameters like Ds,26 will have a smaller dependence on axis 
ratio (V94, Choloniewski 1991). In the same way scale lengths, 
as derived from IRAS-CPC images, can be better understood 
(Van Driel 1995). This way, a randomly oriented galaxy would 
have tb, 21 ~ 1. 

How does this result compare to the current literature? In 
our approach, we compare light which suffers little from ex- 
tinction (K-band) with light which is strongly affected by ex- 
tinction (B-band) within individual galaxies. This is a more 
direct way of deriving the extinction as compared to a sta- 
tistical approach, in which only light affected by extinction is 
investigated and the extinction is derived from studing galaxies 
of different inclinations. We find that at least the galaxies with 
inclinations < 45° are optically thin at Djf,2i or Ds,25- But 
more inclined galaxies could well be opaque at that position, 



so that there is no conflict with Valentijn (1994) or Bottinelli 
et al. (1995), who for large samples in B show that diameters 
at hb = 25 mag (arcsec) - do not vary with inclination. Our 
results are in agreement with Giovanelli et al. (1994) who show 
that at two or three scale lengths spiral galaxies are transparent 
in I. Using the galactic extinction law (see Knapen et al. 1991; 
Jansen et al. 1994) tb/tk ~ 2.2, implying that for our sample 
ti at Djf,2i is at most 0.2, in agreement with Giovanelli et al. . 
They also show that the scale lengths increase with axis ratio, 
more or less in agreement with the present paper. 

We find here that there is in general a considerable amount 
of extinction in spiral galaxies in the central regions, all the way 
out to r e . In the outer parts of our data range, at Djf,2i , or at 3 
stellar scale lengths, we find a maximum optical depth of 0.5 in 
B for a face-on galaxy. Assuming that the dust is distributed 
radially in the same way as the stars, this upper limit would 
go down to 0.1. We find that these galaxies in general, except 
for some edge-ons, are not optically thick at Djf ,21 , although 
some dust could still be present there. 
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